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Abstract 



t^^ . Local duality estimates for the /j-meson leptonic width and the spectrum of radial 

excitations are updated. New perturbative corrections are included in the analysis 



that gives good agreement with experimental data on low energy resonances. 



The method of dispersion sum rules for spectral functions is a powerful tool for study- 



/\ • ing the low energy properties of the hadronic spectrum. The method is based on the 



analyticity properties of the correlators of the interpolating quark/gluon currents and 
asymptotic freedom of QCD. The technique of sum rules was successfully used for de- 
scription of hadronic spectrum in different forms ||I[. Taking into account condensates of 
some local operators within the operator product expansion (OPE) allows one to fix the 
scale of the problem through internal quantities of QCD connected with the properties of 
the vacuum and to predict masses of ground states in a number of channels . 

Within last decade much progress has been made in experimental investigation of low 
energy spectrum of hadrons. Masses of hadronic resonances and their couplings to quark 
currents are now known with better accuracy. The very existence of some states is now 
established more reliably while some of entries have been removed from particle list ^ . 



Theoretically, further corrections in the strong coupling constant expansion have been 
obtained for perturbative parts of current correlators with a variety of quantum numbers 
that calls for updating an old analysis at the higher level of precision. 

While the global duality assumption within OPE is well suited for description of the 
ground state resonances, more detailed information on the spectrum can be obtained by 
assuming the validity of the local duality approximation. The accuracy control is less 
direct in this case. In fact, only an agreement with data justifies the use of the approxi- 
mation, i.e. that the length of the averaging interval is sufficient for comparison between 
quark-gluon and hadronic quantities. Local duality estimates of the characteristics of 
mesonic spectrum happened to be rather successful both for light and heavy quarks 0]. 
Now corrections of higher order in strong coupling expansion for the correlation function 
of the current of massless quarks are available that makes possible to improve on the 
accuracy of theoretical predictions. 

In the present note we consider local duality prediction for the leptonic width of the 
p-meson and masses of its excitations. 

The interpolating current for the p-meson has the form 

which is a neutral component of the isotopic triplet (J^*" = 1 ) and is chosen because the 
resonance residue to the current can be expressed through the leptonic width T{p ^ e^e~) 
measured experimentally. 

The corresponding correlation function reads 

U,,iq) =tJ{Tj,ix)MO))e'^^dx = {q,q, - g,.q')n{q'). 

It is transverse because of conservation of the current. We also consider quarks u and d 
to be massless. For the Adler's function that is free of overall subtraction constant one 
has 

D{Q') = -Q'^n(Q^), -q' = Q\ 



The operator product expansion for D{Q'^) with contribution of the gluon vacuum con- 
densate reads 

D{Q')=D^^{Q^) + {^G^)^^+... (1) 



where 






is the perturbation theory part of the correlator and 

asiQ) 



a{Q) 
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is the couphng constant of strong interaction. Here /ci 2 are known coefficients. 

The function D{Q^) is represented by the dispersion relation with the spectral density 

Pis) 

p{s)ds 



re 

Dm = Q' / 
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The perturbative part of the D-function gives the following expression for the spectral 
density 

+ {h - yPo + i^f^oh + /3i)L + (3^^L'')a%fi^) + . . .) 

with L = ln(yU^/s), /?o,i are coefficients of the /3-function. The term — 7r^/5Qa^(/i^)/3 stems 
from analytical continuation in time-like region of momenta. Such terms can be resummed 
in all orders of a{fi^) for any given order of the /5-function P, P]. 
For fixing normalization we give the definition of the (3 function 



da 

dj? 



^^^n = /5(«) = -/5oa' - Aa' - (32a^ -... (2) 



with 



n ^ n A a 3863 

and numerical values of coefficients ki^2 

fci = — - 9C(3) ^ 1.64, 
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^2 = ^^ - — C(3) + y C(5) - 6.37 

in the MS scheme for Nc = Uf = 3. The coefficient ^2 has been recently computed 
independently in ^ that conffimed the old result [Q]. 

Dispersion sum rules for the spectral density can be based on global or local duality. 
Global duality uses the procedure of smearing for the entire positive semiaxis while the 
assumption of the local duality is that the smearing can be done for parts of this semiaxis. 
Thus the local duality assumption leads to the finite energy sum rules (FESR) of the form 

r ff\s)s''ds = r p'^{s)s^ds (3) 

for some integer k and the integral is performed around the single resonance, p^^{s) 
and p*'^{s) being the physical and theoretical spectral densities correspondingly. Powers 
of energy (s^) can be substituted with more general functions with proper analytical 
behavior. The actual accuracy of the relation (^) depends on the shape of the weight 
functions. Functions without sharp changes are supposed to give the best results and we 
use A; = 0, 1. For intervals not including the origin the set k = —1, is also used. 

Perturbative analysis (without vacuum condensates) of the p meson spectrum within 
the local duality approach was performed in 0] with the result Sq = 2im? and in the 
leading order in as 

2 

T\p^e+e-) = ^-^ = AAkeV (4) 

37r 

where a is the fine structure constant, a = 1/137, and nip = 768.5 ± 0.6 MeV p. This 

result can partly be improved by using perturbative corrections. With strong interaction 

corrections included in the leading order, the old result of is 

Tip^e-e~) = ^l + ^f^^)=4.8keV (5) 
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with A = 100 MeV that was advocated in 0]. The estimate (^ is considerably smaller 
than the present experimental value of the leptonic width |^ 

r^+^^_ =6.77 ±0.32 iceV. (6) 



In the present paper we correct the result (|^) adding the nonperturbative corrections 
(the contribution of gluon condensate) and taking into account new terms of a^ expansion 
for the coefficient function of the unity operator in eq. (|l]). At the level of precision of 
5-10% the running of the electromagnetic fine structure constant (a) that enters the 
definition of the leptonic width eq. (Q) through the residue of the resonance with the 
quark electromagnetic current to the energy scale of the order of 1 GeV is also important. 
All these corrections move the result (|^) in the right direction and give accurate theoretical 
prediction in agreement with experimental data on leptonic width and masses of radial 
excitations. 

Nonperturbative contribution due to nonzero vacuum condensates can be easily taken 
into account 0. For the ground state one has the system of equations for A; = 0, 1 

F = so{l + a{ml)) 



s2 



Fml = |(l + a(mJ))-A (7) 

with the physical spectral density of the form of a narrow resonance approximation 

ff''{s) =F6{s-ml). 
The parameter A 



A=^{^G') = OMGeV' 



r2 
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is a contribution of the gluon vacuum condensate where the standard numerical value 

{^G^) = 0.012 GeV^ is used [|. 

The quantity a is an effective charge for the moments. To make the running coupling 

constant integrable at small momenta we use the technique of effective charges [|l^ . Then 

in the third order of PT for the /3 function there is an infrared fixed point JlTI. The results 

of integration are 

/ dsa{s) = aoSo, / dssa{s) = ai — 

Jo Jo 2 

with ao = 0.302, ai = 0.281 for sq ^ 1.3 GeV^ and and a{ml) = 0.12(0.380/7r) [|l3 



chosen as an input, tti,- = 1.777 GeV is the r lepton mass. Note that the semileptonic 



width of the r lepton that is now one of the best sources for low energy determination of 



the couphng constant is given by a particular FESR [T^. In the following we neglect the 
difference between oq and ai putting oq ~ ai ~ a = 0.3. 

Because both A and a are small we treat them as corrections and limit ourselves to 
the linear approximation in these parameters. 

From eqs. (|^) one has 

2 A 

So = 2m H 

^ rrip 

and 

T^+^^ =r°(p^e+e-)(l + a)(l + -) = 4.4 * 1.3 * 1.06 = 6.06 ieV 

2m p 

that is still smaller than the experimental number (^). First factor stems from perturbative 
corrections and second from gluonic condensate. 

An additional contribution comes from remormalization of the electromagnetic cou- 



phng constant (for details, see e.g. WW) 
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1 — Aa(s) 37r ^ ■' y m^ 



f 
where the sum runs over the fermions with masses smaller than s, Qj is the fermion 

electric charge, N^f = 3 for quarks and Ncf = 1 for leptons. We find 



1 - Aa{1.3 GeV^) 1-0.025 
that gives another five percent in the width. Finally we have 

r(p^e+e-) = 6.37]feV (8) 

which is the main result of the paper that is much closer to the experimental value (P) 
though is still somewhat low. The agreement can be achieved by increasing the numerical 
value of the gluon condensate that is not very plausible because the masses of excitations 
will be larger. 



The empirical estimate 



ml + ml, „ , ^ 

So = ^ ^ ^ = 1.35 GeV^ (9) 



also gives a good agreement with the experimental value, namely 

T^^p(p _ e+e-) = 6.76 keV (10) 

because the duality interval (H) is somewhat larger. 

Next, the radial excitations can be predicted then using the assumption that the 
integration border lies exactly in the middle between subsequent resonances |Q. 

Corresponding parameters for p' are ml, = 1.96 GeV^ {nipi = 1.40 GeV) while the 
experimental number is 

nip, = 1449 ± 8 MeV 

with the total width T = 310 ± 60 MeV. 

It is reasonably good numerically though is not very well justified just because of 
large width of the corresponding state. The point is that the very definition of mass of 
the resonance depends on the procedure used for fitting data. For instance, there is a 
difference whether the common Breit-Wigner resonance curve is used in terms of a center 
mass energy E or the energy square s = E"^ (see a thorough discussion in ||14|| ). 

Even higher excitation is in a reasonable agreement with our naive formula, namely 

^P(i700) = 1717 ±13 MeV 

while we predict 

mp(3) = 1.80 GeV. 

For the next one we have mp(4) = 2.13 GeV while the experimental number is 

mp(2i50) =2149 ±17 MeV 

These local duality estimates of the spectrum that are essentially the linear model of 
ref. m are in better agreement with experiment because experimental situation changed. 
The state p(1250) moved top (1450) andp(1600) is shifted to p(1700)(m = 1717±13Mey) 
improving agreement with local duality predictions. Still for the p(1700) the discrepancy 
is rather large. 



Note that for exited states one can compute the resonance parameters using two sets 
of sum rules for k = 0,1 and k = —1,0. The difference between predictions can be 
considered as a practical estimate of the accuracy of the method. Roughly, numerical 
estimate is needed in any particular case but for two sets mentioned above the difference 
is within 15%. 

Main ingredients of the improvement of the agreement of finite energy sum rules 
predictions with experimental data are: 

• Better experimental data, especially for first two excitations. The state p(1250) 
moved to p(1450) that allowed to increase the duality interval for the p- meson 
without contradiction with the spectrum. 

• Better knowledge of the QCD coupling and understanding of its infrared behavior 
that gives the agreement for the width but practically does not affect the ratio of 
the moments for determination of the duality interval and the mass of the first 
excitation 

• An account for the gluonic condensate that moves numerical estimates in the right 
direction. Its effect is rather small for the standard numerical value of the conden- 
sate. However it improves on the agreement. Note that it cannot be increased much 
if one would like to preserve the spectrum of radial excitations within this approach: 
for larger condensate width fits better the experimental value but predictions for 
masses move too far. One should not be too strict on that however because for 
higher states the approximation of narrow resonances may become inadequate and 
finite width of resonances should be taken into account. 

• Renormalization of the fine structure constant also improves the prediction for the 
p-meson leptonic width. 

Note that perturbative corrections cancel in the ratio for determination of the duality 
interval through the mass of the ground state. Thus the prediction for the mass of radial 
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excitations is not sensitive to pQCD corrections while the residue is directly proportional 
to their magnitude. Their account improves the agreement with experimental data. 

Agreement with experiment is rather good. That makes FESR a powerful tool for 
working out the phenomenology though the internal mechanism for determination of the 
ground state mass is absent and the mass of the ground state sets the scale in deter- 
mination of the spectrum. Thus these two approaches are complementary. As far as 
condensates are used to find a mass of the ground state the FESR can be used for further 
investigation of the spectrum. Power corrections do not affect higher radial excitations 
within the local duality approach. 

The low energy states are now much better known that allows to reanalyzed some 
FESR results with other quantum numbers as well. 
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